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Abstract 

By a semi-empirical Green's function method we calculate con- 
ductances and the corresponding Giant Magneto-Resistance effects 
(GMR) of two metallic ferromagnetic films separated by different spac- 
ers, metallic and non-metallic ones, in a simplified model on a sc lattice, 
in CPP and CIP geometries (i.e. current perpendicular or parallel to 
the planes), without impurities, or with interface- or bulk impurities. 
The electronic structure of the systems is approximated by two hy- 
bridized orbitals per atom, to mimic s-bands and d-bands and their 
hybridization. 

We show that such calculations usually give rough estimates only, 
but of the correct order of magnitude; in particular, the predictions on 
the impurity effects depend strongly on the model parameters. One 
of our main results is the prediction of huge CPP-GMR effects for 
non-metallic spacers in the ballistic limit. 
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1 Introduction 



The presence of giant magneto-resistance effects in trilayer systems 
consisting of two ferromagnetic metallic layers - source and drain - 
separated by a non-magnetic or anti-ferromagnetic crystalline metallic 
spacer is meanwhile well-known [p], |2]. Recently also non-crystalline 
and non-metallic spacers have been considered ||, |], |5[, i.e. there is 
renewed interest in spin-polarized tunneling through semi-conducting 
or non-conducting spacers || [7], |§, based on the perspective of new 
applications in magneto-electronic technology, e.g. for magnetic-field 
sensors, spin-valve transistors || or spin-polarized field-effect transis- 



tors \m li 



In |12|, 13 spin- valve properties of ferromagnet /insulator /ferromagnet 



junctions were studied, and presents a first ab initio calculation 



of the one-particle electronic and magnetic properties of FM/NM/FM 
tunneling structures, with FM=Fe, and NM=Ge and GaAs. 

However concerning the GMR effect, which is based on two-particle 
properties, equally accurate calculations for FM/NM/FM systems, par- 
ticularly in the CPP geometry, apparently could only be performed 
without impurities, |T5| . For systems with metallic spacers one should 
mention at the same time the calculations of Tsymbal and Pettifor, 
|TE| |, or Mathon, |T7|] , where also fully realistic calculations have been 
performed, but again with phenomenological assumptions for the im- 
purity scattering, [HJ, or with no impurities at all, [17]]. Furthermore, 
in a two-band tight-binding approximation, two of the present authors, 
have already treated systems, where one or two of the ferromagnet- 



spacer interfaces were decorated with ultrathin non-metallic layers, but 
again with translational invariance within the planes. In this approxi- 
mation, the two bands mentioned above describe the s- and d-electrons 
and their hybridization, or the conduction and valence band in case of 
a non-metallic spacer. 

In the ferromagnetic metals, the parameters of the d-band are of 
course spin-dependent. 

Since in [18] the method for the calculation of the resistivities turned 
out to be extremely flexible and accurate, we extend it in the present 
paper to a systematic survey of the magnetoresistivity of trilayers of two 
ferromagnetic sandwiches separated by different non-magnetic spacers, 
metallic and non-metallic ones, with impurities of various kind, and 



2 



without impurities, in CIP and and CPP geometries, but within a 
simplified semi-empirical hybridized two-band model. 

We stress that our approach is rigorous in principle, numerically 
very accurate, and avoids the Coherent Potential Approximation (CPA), 
which is often applied to disordered systems, but difficult to control. 

As already mentioned, there exists already a large body of theo- 
retical calculations on the GMR of magnetic multilayer systems. For 
systems with metallic spacers and without impurities our simplified 
two-band model calculations, to be presented below, can of course not 
compete - and this is also not at all intended - with the much more 
ambitious and cumbersome calculations mentioned above e.g. those of 



17 1, who performed fully realistic calculations of the CPP- and CIP- 



GMR for Co/Cu/Co/Cu and Fe/Cr/Fe/Cr multilayers with ideal leads 
made of Cu rsp. Cr, in the ballistic, Ohmic, and Anderson-localized 
regimes corresponding to the cases of no-disorder rsp. certain cases of 
one-dimensional disorder, i.e. with respect to the stacking of the planes, 
but without impurities. There exist more such ambitious calculations, 



e.g. the seminal paper of Schep et al, [[L9), and other publications, 
which also in principle would have earned citing but which are omitted 
on lack of space. In this respect the purpose of our paper is, instead, 
to study in some detail the question, to which extent simplified models 
as the present one are able (or not) to reproduce those results. 

On the other hand, with impurities, on which we concentrate, sim- 
ilar extensive calculations of corresponding rigour do not yet exist, to 
our knowledge. In fact, concerning the influence of impurities, the 
methods of our calculation have already been presented in the paper of 
Asano et al, p0| . However, rigorous calculations with impurities are so 
demanding that those authors, although they wrote already down the 
formalism for a two-band model, published numerical results only for 
one-band cases. 

Here we present results obtained with the two-band model within 
the formalism of f20fl . For a nonmetallic spacer it is actually neces- 
sary to consider at least two orbitals per atom to simulate the valence 
and conduction bands, respectively, but also for ferromagnetic or non- 
magnetic 3d-metals it is a natural requirement to include a second band 
with the intention to mimic as far as possible the d-electrons in addi- 
tion to the s-states. Namely, the d-electrons are not only responsible 
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for the magnetic properties, but also for a large part of the resistivity, 
due to Mott's s-d-scattering mechanism. Therefore it is also impor- 
tant to include not only two bands, but also realistic values for the s-d 
hybridization. 

This should give enough motivation to our two-band calculation; 
but we stress again that these calculations should be considered as 
model calculations and are not intended to replace fully realistic calcu- 
lations, as far as they are possible already now, e.g. |T7j, or in future. In 
fact, it turns out that the present two-band model calculations will be 
able to give estimates of the GMR, which yield the correct order of mag- 
nitude, but not more. Further, concerning the impurity effects, it turns 
out that certain contradictory results (see below) concerning essential 
trends depend strongly on the model-parameters, so that again, now 
with impurities, completely realistic calculations cannot be avoided: 
Unfortunately, such calculations, where the impurities are treated with 
the same kind of rigour as the electronic structure of the pure system 
itself, are apparently not yet possible, f21"|j-f25|j. 

On the other hand, one of our most remarkable predictions, which 
should not be overlooked and can be stated already from the present 
model calculations, is that a large enhancement of the CPP-GMR may 
happen, if nonmetallic spacers of increasing thickness are used. Here 
it is of course required that the thickness of the nonmetallic spacer 
remains smaller than the dephasing length for inelastic scattering; i.e. 
one should remain in the ballistic limit. 

In the following sections we present the basic theory and then the 
results of our survey, which is followed at the end by a section present- 
ing our conclusions. 



2 Basic definitions 

2.1 Systems 

Our systems are defined in Fig. 1 and Fig. 2 for the CPP and CIP 
geometries, respectively. 

We start with metallic 'ideal leads' on the left-hand side and right- 
hand-side, respectively, followed - in the CPP geometry - by three 
monolayers of ferromagnetic metal F% rsp. F%, with n s spacer mono- 
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layers inbetween. In the CPP geometry all layers have a quadratic 
cross-section of M x M square-lattice unit cells. 

For the CIP geometry, the natural definition of our systems is sim- 
ilar and explained by Fig. 2. 

The electronic structure of our systems is described by the following 
tight-binding Hamiltonian : 

ft = Y E ti Y ^tAi.cr + Y Y tS l,i,V,i' Y CM,<T^',j',<7 
Li <? Li V A! cr 



+ Y Y E l,i,a d%i,trdl,i,<r + YY *M,i',*' Y ^i^'.iV 
l,i & l,i V,i' CT 

+ Y V lf Y^t^i,* + Cl M dl,i,a) . 



li <? 

Here o[ i a and denote creation and destruction operators of 
an s-electron occupying the site I in plane i with spin a. The corre- 
sponding operators for d-electrons are a and d\^ a . The band-energy 
parameters are Ef i for the s-states, which do not depend on the spin 
a, and Ef i(T , which are spin-dependent, for the d-states. The hopping- 
matrix elements for s- and d-states are tf i v v and tfj v v , and the local 
s-d-hybridization is given by V t s - d . 

Of course, taking more bands into account, i.e. all five d-bands, the 
4s-band and the three 4p-bands, and adding the spin dependence, we 
could have tried to fit the band structure of realistic systems as far as 
possible, see e.g. [plf] ; but this is notour purpose, since we try to remain 
semi-quantitative and study only the essential trends. Particularly, in 
view of the still much simplified band structure, we also use typical 
simplifications in the description of the impurities, see below. But in 
principle, within our formalism, more rigour in both respects would be 
possible, but only at the cost of excessive computing. 

For the following, H. is written in matrix representation with respect 
to an orthonormal basis by replacing the creation and annihilation 
operators by the corresponding ket- and bra-vectors, respectively, such 
that TC — » H , e.g. by d^ a h,io \d, /, i, a) (s, I, i, a\. 

2.2 Resolvent operators and Kubo formula 

For the calculation of the conductances we use the formalism of Fisher 
and Lee, |27[|. Therefore, we need the following matrix elements of the 
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advanced and retarded resolvent operators G ± := {E ± i0 + — H} 1 : 
G^M'WjoW = (a,i,l,a\{E ±i0+ -H}- l \a',i',l', a') . (1) 

Here the index a counts the orbitals, e.g. a = s or d. As usual, we also 
call these matrices 'Green's operators', and i0 + is a positive imaginary 
infinitesimal, such that the matrix-inverse in this equation always exists 
in our approach. 

For auxiliary purposes we need the 'left-sided' and 'right-sided' 
Green's operators G L and G R , e.g. 

G L (t ) a ,i,a;a>,i>,a> ■= (a, i , 1, a\{E + i0 + - H?}' 1 ^' , i , V , a') . (2) 



Here denotes the Hamiltonian for a system, in which all planes 
i > iq are deleted; G R (io) is defined in a similar way with E+i + — , 
where is the Hamiltonian of a system with all planes % < %q deleted. 
G L and G R are obtained by recursion, e.g. G L by recursion from the 
left, i.e. 

G L (i ) = \g(io)- 1 -fG L {i -l)f\- 1 . (3) 

g(io) is the resolvent of the isolated plane io] T is the matrix describing 
the 'hopping' from plane to plane. 

From this one gets finally the two desired resolvent operators 

G(i, i) = [git)- 1 - TG L (i - l)f - f G R (t + l)f ]- x (4) 

and 

G(i,i + l)=G(i,i)TG R (i + l), (5) 

which are needed for the conductance, see below. Here the plane % 
(= io) is arbitrary, since under stationary conditions through any of 
the planes the same amount of current is flowing. 
Finally, with the definition 

G(i,i') : =I[G-(z,0-G + ( ? ,0] (6) 



one gets the conductance from the Kubo-formula, [g? 



T a = -^Tr, a [G{i,i)fG(i + l,i + l)f-fG(i,i-l)TG{i,i-l)}, (7) 
h 
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where the trace Tr^ Q is performed with respect to the orbitals a = s, d, 
and with respect to the sites I of the atoms of a given plane, whereas 
the spin-index o (= a') is kept fixed. Here e is the electronic charge, 
and h is Planck's constant. 

The GMR effect is defined by the equation 

y++ i p++ -p++ 

GMR = -2=1 -pi - 1 = - 1 . (8) 

pH — I pH — pn — v ^ 

Here the superscripts (++) resp. (H — ) denote mutually parallel rsp. an- 
tiparallel magnetizations of the two ferromagnetic sandwiches. Since we 
neglect spin-flip scattering altogether, the conductances are of course 
additively composed out of the separate contributions from majority 
and minority spins, respectively, as denoted in the formula. 

We have performed real-space (x-space) and k-space calculations. 
In the x-space calculations each layer consists of M x M =10x10 
atoms with free boundary conditions, PBfl . Here the matrices are of 



size 200x200, again for given a (= a'). Impurities can be added at 
will, which is important below. 

In our A;-space calculations, which are only performed in the CPP 
case and without impurities, the planes are infinitely extended. Here, 
in the trace of equation ([/]), in the fc-space representation, the sum 
over / is replaced by the corresponding sum over k- vectors in the two- 
dimensional Brillouin zone, since the two-dimensional vector ku is now 
'a good quantum number', i.e. there is translational invariance in the 
planes. This summation, or the corresponding integral, is approxi- 
mated extremely accurately by a Cunningham formula, [p9|1 , with ~ 10 6 



fc||-points, for which the matrix elements of the Green's operators can 
be calculated separately without difficulty and extremely accurately. 
(All numerical evaluations have been performed with Cunningham's 
accuracy parameter m = 10 or 11). 



2.3 Model parameters 

For the models considered, we have chosen four parameter sets: Two 
parameter sets, (i) and (ii), correspond to a metallic spacer and mimic 
the main features of trilayers composed of Co as ferromagnetic metals 
and Cu as nonmagnetic metallic spacer material (see below), whereas 
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the third and fourth parameter sets, (iii) and (iv), correspond again to 
Co ferromagnets, but with a non-metallic spacer. 

Concerning the sets (i) and (ii), we note that, as already mentioned, 
they are constructed to reflect the two main features of Co and Cu as 
far as possible in our simplified two-band model: The two features are 

• the similarity of the s-bands of Co and Cu, and 

• the similarity of the d-bands of Cu and the 'majority-spin' d- 
bands of Co. 

This qualitative similarity is obvious from Fig. 3, which presents results 



obtained by Mathon et al, |30| . The parameter sets (i) and (ii) are 
in fact derived from rough fits to the overall band-structure of Co - 
parameter set (i), see Fig. 4a - respectively to the overall structure of 
the DOS for Co - parameter set (ii), see Fig. 4b. But these fits are 
rather crude anyway, and so the derivation is not emphasized. The 
values are given in Table 1 and Table 2, respectively, [ |31|j . 

In Tables 3 and 4, the parameter sets (iii) and (iv) are given, which 
correspond to nonmetallic spacers ('Isolator 1 and 2', respectively). 
The parameters of the ferromagnet correspond to case (i), with Ep = 
-2.8. 

The notions 'indirect energy gap' and 'direct energy gap', which are 
used in the table captions of Table 3 and 4, respectively, refer to Fig. 5: 
For parameter set (iii) there is an 'indirect energy gap\ since the valence- 
band maximum is at fc[m] = (ir/a)y3 with an energy -2.81, which is 
slightly below the Fermi energy Ep = —2.80, whereas the minimum 
of the conduction-band energy is at fyni] = with E=-2.79. These 
numbers refer to the first line of the three parameter sets of Table 3. 

On the other hand, with parameter set (iv), also with Ep = —2.8, 
there is a direct energy gap at fcmii = 0.814, i.e. at the arrow in the 
r.h.s. of Fig. 5. In both figures, and in the results presented below, the 
uppermost line, i.e. with the smallest gap, of Table 3 and Table 4 has 
been used. 

2.4 Impurity models 

We consider a) bulk impurities and b) interface impurities. 
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In case a), we have considered bulk impurities only in the spacer, 
and exclusively nonmagnetic ones, which are modelled by adding to 
the parameter E d of the single-site d-band energy a Gaussian spatial 
noise proportional to a 'disorder strength' oy, namely 

E d - E d + a r ■ np auss . (9) 

Here the np auss , where I enumerates the different spacer atoms, are 
independent random numbers distributed according to a 'Gaussian' 
with average and variance 1; the parameters of the s-bands remain 
unchanged. 

These assumptions for our bulk impurities are of course rather 
schematic, i.e. these are the same simplified assumptions e.g. of Asano 



et ai, [20]. But principally one might have been more ambitious and 



might have considered e.g. specific magnetic impurities. Such calcula- 



tions are of course always possible within our formalism, [32. 33], and 
have partially been performed in course of the diploma work of the 
first one of the present authors; however, as a first step, and since the 
results with the magnetic impurities were not particularly exciting and 
have not been documented, also because of lack of space, we keep at 
present to our simplified description of the impurities. This seems also 
justified in view of the simplifications made for the band structure and 
for the geometry of our systems. 

Concerning case b), the interface impurities are defined as follows : 
For every ferromagnet-spacer interface, in each of the two boundary 
layers of the interface there are n atoms (out of M x M) selected at 
random and replaced by atoms of the other kind. 

For impure systems we have of course averaged over a large number 
of different samples, see the figures below; this leads to the 'averages 
with error bars' in the results presented in the following section. 



3 Results 

3.1 CPP-GMR 

3.1.1 CPP-GMR; no impurities 

At first we describe the results obtained without impurities obtained 
in the CPP case ('current perpendicular to the planes'). In Fig. 6 
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and Fig. 7 the results of k-space calculations with model parameters 
(i) and (ii) are presented, for variable thickness n s of the 'metallic 
spacer'. Interestingly, in both cases there is a slight but significant 
spatial oscillation of the GMR as a function of n s ; in case (i) the GMR 
oscillates between ~ 70 and ~ 80 % with a quasiperiod of A ~ 5 
monolayers (ML); in case (ii) the oscillation of the GMR is between 
~ 47 and ~ 51 % with A ~ 3. Here one can explicitly see that the 
conductance is ballistic from the fact that on average the conductances 
hardly depend on n s in the limit n s — > oo; in contrast, in case of 
Ohmic behaviour, one would expect that the conductances converge 
to 0, namely oc nj 1 , as soon as n s becomes larger than the relevant 



"scattering length" for diffusive elastic or inelastic scattering, |17 . 

The reason for the spatially oscillating behaviour of the GMR are 
tiny oscillations of the conductances r+ + and r^ + , which are visible in 
the plots on the r.h.s. of the preceding figures. These oscillations are 
caused by our treatment of the ideal leads: In fact, for the ideal leads 
we have always used the same parameters as within the 'Cu' spacer, 
but usually - with exceptions mentioned below - for the pure system 
the s-d-hybridization has been neglected in the ideal leads (but not in 
the spacer). The reason for this neglection is mainly technical, since 
in case of neglected hybridization we can directly express the matrix 
elements of the Green's operators G L (k, i) and G R (k, i) at the boundary 
layers from the ideal lead to the ferromagnetic by simple analytical 
expressions, whereas the numerical calculation is much harder if the s- 
d-hybrization is taken into account already in the ideal leads, by which 
the above-mentioned tiny oscillations vanish. What is, however, more 
important: In this way we have made explicit that the properties of 
the ideal lead really matter, which is not astonishing with 'ballistic 
electrons', where the 'contact resistance' from the ideal leads to the 
system under study plays an important role. 

As a first consequence one can learn from this study that it is a mis- 
take to speak of 'the ideal lead' as a uniquely defined entity: In ballistic 
experiments there are different ideal leads, and it is necessary to char- 
acterize them, too, as already stressed in a recent book on electronic 
transport in mesoscopic systems, |34[ . 



It is also interesting at this place to contrast the different, i.e. 
parameter-dependent, results of our Fig. 6 and Fig. 7 with the cor- 
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responding results found by the much more ambitious ab initio calcu- 
lations of [|T7| and ||19|| . These authors find values of the CPP-GMR 
for Co/Cu- multilayers amounting to ~ 150 %. Moreover, for the con- 
ductances T/M 2 [e 2 /h], instead of our values 0.8; 0.3; and 0.2 (Fig. 6) 
rsp. 0.45; 0.2; and 0.15 (Fig. 7) for the cases (++,T), (++,!), and (+-), 
respectively, they find values of T corresponding to 0.43-10 15 , 0.25-10 15 , 
and 0.165-10 15 Q" 1 m~ 2 , i.e. T/M 2 [e 2 /h] = 0.67, 0.39, and 0.26. 

As a conclusion, our values are correct concerning the order of the 
three conductances, and their order of magnitude, but not more. This 
is exactly what can be expected. 



3.1.2 CPP-GMR; bulk impurities 

We would like to stress at this place that with impurities, all calcula- 
tions have been performed in the real lattice space, with cross sections of 
MxM =10 x 10 sites, free boundary conditions, in the perpendicular di- 
rections, instead of periodic ones, and with complete s-d-hybridization 
everywhere, i.e. now always including the ideal leads. 

In Fig. 8 and Fig. 9, for n s = 3 and the parameters (i) and (ii), re- 
spectively, our results for the GMR and the conductances are presented 
over the strength of the disorder, i.e. over the standard deviation ay 
in equation (|9]). For the parameter set (i), Fig. 8, one obtains first a 
rapid decrease of the GMR from GMR « 1.15 to GMR ps 0.45, when 
o~ r increases from to 1, followed by a plateau behaviour with further 
increase of ay. Concerning the behaviour of r+ + , r~~, and r+~, one 
can see from the r.h.s. of Fig. 8 that the GMR reflects essentially the 
behaviour of 

For the parameter set (ii), Fig. 9, there is also at first a rapid de- 
crease from GMR « 0.37 at o r = down to GMR « 0.15 at a r » 1, 
but with a further increase of a r there is now a significant re-increase of 
GMR back to a r « 0.22 for o~ r « 6. In fact, in this case the behaviour 
of the conductances, which are plotted on the r.h.s. of Fig. 9, is sig- 
nificantly different and more subtle than that of Fig. 8. In particular, 
the results of Fig. 9 differ from what in |2(J the authors have obtained 
from a simple one-band model. 

So from these two figures, Fig. 8 and Fig. 9, concerning the de- 
pendence on the strength of disorder o r , one may only conclude that 
at first a rapid decrease of GMR(a r ) from its value at a r = to 
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GMR(l) pa 0.4 • GMR(O) would be typical; this qualitative behaviour 
seems to be essentially parameter-independent; but for a r > 1 a quan- 
titative calculation is hardly avoidable, |22|, |55| . 

The reason for the common 'decreasing behaviour 1 for a r <l is the 
following: In this case, an increase of o~ r leads to an effective reduction 
of the number of open channels for coherent motion of an electron 
from one plane to the next. This reduction concerns mainly T+ + and 
therefore leads to a decrease of the GMR, since the majority-spin d- 
electrons, without impurities, would not feel any potential- variation at 
all. 

In this connection, a very recent paper of Bruno et ai, ||36|| , should 
be mentioned: There it is shown that the conductance of systems with 
impurities can be decomposed into two parts, a 'diffusive 1 part, aris- 
ing from scattering processes, where the in-plane /c||-vector is not con- 
served, as is generically always the case for the present systems with 
impurities and free-boundary conditions, and a ' 'ballistic'' part, to which 
only those processes contribute, where the k\\- vector is the same for the 
incoming and outgoing waves. The diffusive part ( rsp. ballistic part) 
always dominates the conductance, as long as the spacer length is much 
shorter (rsp. much longer) than the scattering length l s of the system. 



According to |$6] , l s should be of the order of several 100 lattice spac- 
ings under similar conditions as in the present work. So in our systems 
with impurities, where n s is very small (n s = 3), always the diffusive 
part of the conductance dominates. 



3.1.3 CPP-GMR; interface impurities 

In the following Fig. 10 and Fig. 11 the number n of impurities is varied 
between n = and n = 80 for every interface layer (remember that 
per interface, two layers are involved, with 100 sites per layer; so we 
have 2 x n impurity sites out of 200 per interface, occupied by a 'wrong 
atom'). 

For the parameter set (i), Fig. 10, the GMR decreases with increas- 
ing number of impurities, whereas for parameter set (ii), Fig. 11, there 
is again the opposite behaviour: Here the GMR increases with increas- 
ing number of impurities. 

This difference contains again an important statement in itself, which 
should be particularly relevant to the experimentalist : 
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Without a precise calculation for a specific model one cannot predict 
whether an increase of the number of interface impurities will lead to 
an increase or a decrease of the GMR. 



3.2 CIP-GMR 

3.2.1 CIP-GMR; no impurities 

The formalisms of the preceding sections can also be applied on equal 
footings to the situation, where the current is in the direction of the 
planes (CIP geometry, see Fig. 2). The ideal leads have cross-sections of 
AzxAy = NxM with N = 6+n s , see below, while the CIP-conductor 
consists again of Az = N monolayers of length L = Ax = 10 and width 
Ay — M = 10, namely two ferromagnetic sandwiches of heigth Az = 3 
and inbetween a spacer sandwich of Az = n s . Again we consider at 
first the situation without impurities; however in this geometry we have 
applied the x-space calculation right from the beginning, i.e. with cross 
sections of size 10 x 10 and with free boundaries, although with respect 
to the coordinate y one could have still performed a Fourier transform. 
Moreover, as for the CPP geometry without impurities, this time the 
s-d-hybridization has been switched off for the lead wires. 

In Fig. 12 and Fig. 13, the ballistic CIP conductances are plotted as 
a function of n s , for the parameter sets (i) and (ii), respectively. In the 
first case, with increasing n s the conductances T/M 2 rise, because the 
spacer cross section increases, the limiting value given by the "shunt- 
ing" through the Cu spacer. The GMR shows irregular oscillations as a 
function of n s ; these oscillations change with the boundary conditions, 
the geometry, and the material, and the results are more than one order 
of magnitude smaller than for the CPP geometry, |3S|, |3I|. Note 



that with parameter set (i) , the GMR is even negative on average over 
n s , which is not observed experimentally. By discrete Fourier transfor- 
mation of GMR with respect to n s , [flOfl , we would get maxima in the 
Fourier spectrum at periods of 2.5, 4.6 and 8.33 monolayers (ML) from 
parameter set (i), Fig. 12, and 2.27 and 5 ML from (ii), Fig. 13. 

3.2.2 CIP-GMR; bulk impurities 

In case of bulk impurities, we have considered a system with n s = 4, i.e. 
iV = L = 10. The results are presented as a function of a r , the param- 
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eter characterizing the 'strength of the disorder' according to Eq. (£5|). 
Here, in case of the parameter sets (i) and (ii), we find again quite dif- 
ferent behaviour in Fig. 14 and 15, respectively: In case of Fig. 14, i.e. 
for parameter set (i), the GMR increases at first from to 0.015, and 
then for o~ r > 1 it roughly remains constant with increasing a r , whereas 
in case of the parameter set (ii), i.e. in Fig. 15, there is at first a drastic 
decrease from GMR ~ 0.035 to GMR ~ 0, and then again a roughly 



constant value of GMR ~ 0, p7B - In both cases the conductances de- 
crease at first rather fast, and then there is a slow re-increase; but this 
behaviour is slightly different for the up- and down-spin channels in 
the (++)-configuration, and for the (H — )-configuration. 

So again the behaviour of the GMR as a function of o~ r is hardly 
predictable without an extensive and completely realistic calculation for 
a specific model. 



3.2.3 CIP-GMR; interface impurities 

A somewhat different conclusion refers to the influence of interface 
impurities for the CIP geometry. Here in Fig. 16 and 17 we use the same 
geometry as in the two preceding figures. The GMR rises significantly 
up to 50 % impurities - i.e. this time we have the same trend for both 
models (i) and (ii) -, and then it falls again, which is natural since 
a system of two three-ML ferromagnets separated by a 4-ML-spacer 
and x > 50 % interface impurities is equivalent with two four-ML 
ferromagnets separated by a 2-ML spacer with x' = (1 — x) < 50 % 
interface impurities. Here it should be noted that the majority-spin 
carriers are hardly influenced by our interface disorder, which is clear 
since we have assumed that the s-bands and the majority d-bands of 
our metallic ferromagnets and of the metallic spacers are identical. 
But it should also be noted that at the above-mentioned maximum the 
CIP-GMR with interface scattering reaches much higher values (e.g. 
9% and 8% respectively, for the parameter sets of (i) and (ii)) as with 
our nonmagnetic impurities in the bulk. 
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4 Non- metallic spacers, CPP geometry 



Finally, we come to our most startling result, for the non- metallic spac- 
ers, where we have considered metallic ferromagnets as in model (i), 
but semi- conducting spacer layers as in the first line of Tables 3 and 
4. The geometrical situation corresponds to the CPP geometry of Fig. 
1, with metallic ferromagnets (two-times 3 monolayers) and a variable 
number n s of semi-conducting spacer layers, i.e. with a conduction band 
and a valence band, inbetween. Of course one now expects a strong 
exponential decrease of all conductances with increasing n s , but GMR 
= (r ++ /T + ~) — 1 may be well defined and even increase with increas- 
ing n s , as long as the spacer thickness does not become larger then the 
'dephasing scattering length' or 'spin-flip scattering length'. 

The calculations have been performed by the highly accurate k- 
space method in reliable numerical accuracy, and we have plotted the 
conductances and the GMR against the number n s of spacer layers, 
calculated in the ballistic limit of perfect phase coherence. In Fig. 
18 and Fig. 19 there is now the expected drastic exponential decrease 
of the conductance, which for our model is weaker for T ++ than for 
(The decay is weakest for r+ + , which is reasonable, since the j 
d-electrons have a much lower energy gap to the lower edge of the semi- 
conductor conduction band than the f d-electrons, and with the T + ~ 
case, in contrast to the T ++ case, there are the additional 'matching' 
problems at the interface to the second ferromagnet. Note that due to 
the s-d-hybridization, and with the finite value of td, the d-electrons 
do in any case contribute to the tunneling current.) So this means 



that in our case, and in the recent calculation of MacLaren et al, |L5 



to be discussed below, and in contrast to Julliere's assumptions, @, 
the tunneling matrix through the semi-conducting spacer is naturally 
strongly spin-dependent. (The fact that here for T" 1 the exponential 
decay with the spacer thickness is stronger than for r ++ , depends of 
course on our model: for other assumptions on the semi-conductor it 
might be just the other- way round.) 

As a consequence, one gets a positive GMR, which strongly increases 
withn s , e.g. up to GMR ~ 300 % for n s > 5 in case of model (iv), Fig. 18, 
i.e. for the model with the 'direct energy gap'. For the case of model 
(iii) with the 'indirect gap', Fig. 19, the increase is at first similar, 
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i.e. for n s <4, but then the GMR increases even further to collosal 
values, e.g. GMR ~ 2500 % for n s = 7. However because of the strong 
exponential decrease of the conductance, and in view of the fact that 
we have neglected indirect transitions involving phonons, this region 
is probably beyond experimental realizability, and the results of Fig. 
18 and Fig. 19 should not be taken literally, although the numerical 
accuracy of our calculations seems to be sufficient (In Fig. 18 and Fig. 
19, the Cunningham accuracy parameter, see |29|], was m = 10). 

In fact, however, in the already cited recent paper of MacLaren et 



ai, | 15| , which contains an ambitious realistic LKKR calculation of the 



ballistic CPP conductances and the corresponding GMR for bcc (001) 
Fe/ZnSe/Fe trilayers, the conductances of Fig. 3 of that paper resemble 
very much - almost quantitatively - to the results of our simple model 
calculation in Fig. 19. One only must take into account that in the 
paper of MacLaren et ai, one Zn-Se bcc double-layer is replaced by 
two monolayers in our simple sc two-band model calculation with the 
parameter set (iii). So here we have an example where a fully realistic 
ab-initio calculation may be sometimes not necessary after all. 

Finally we present at this place additional comments on the sem- 
inal paper of Julliere, 0, written already in 1975. Julliere describes 
the tunneling conductance from the ferromagnetic metallic source to 
a metallic ferromagnetic drain through an isolating Al 2 03 spacer as a 
product of the densities of states (DOS) of the ferromagnetic source 
and drain, multiplied by an effective tunneling matrix element. This 
should be essentially equivalent to our approach, which involves two 
imaginary parts of Green's function, which represent in principle den- 
sities of states. But from the relevant densities of states of the magnetic 
metals, the contributions of confined electrons, which contribute to the 
DOS, but not to the conductance, should be excluded; such confined 
states exist e.g. with the Co down-spin electrons, which are strongly 
reflected at the Cu interfaces; so instead of Juliere's DOS, an 'effec- 
tive DOS' should be used, and since the reduction factor of the DOS 
to this 'effective DOS' is spin-dependent, the remaining 'effective tun- 
neling matrix elements', if one uses the original DOS, i.e. Julliere's 
formula, would in any case depend on the spin, contrary to Julliere's 



assumption, if applied to systems as C0/AI2O3/C0, j|T 



In principle, however, our approach is even more explicit, and more 
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demanding, since it also involves the transfers to (and from) the ideal 
leads to the ferromagnetic metal, and from the reservoirs to the ideal 
leads. 

In Julliere's TMR experiments, and similar experiments performed 
at present, the high values of the TMR predicted in the present pa- 
per and also in the paper of MacLaren et ai, fll5| , have not been ob- 
served hitherto. However one should stress that our insulating layers 
are extremely thin (from 1 to 15 monolayers only) with ideal inter- 
faces, whereas the experimental thicknesses of the AI2O3 spacers are 
larger and with rougher interfaces, so that in that case spin-channel 
mixing might play a role, particularly if on the metallic side of the 
interface the transverse spin components are almost as probable as the 
z-components. The band-structure, e.g. that of the semi-conductors 
involved, is also crucial, since the exponential decay with the spacer- 
thickness should be weaker for T ++ than for r + ~, i.e. the semi-conducting 
compounds involved should have properties similar to those of our 
present models, or to the Fe/ZnSe/Fe-system studied by MacLaren et 



ai, [15|. In any case, one can only speculate that perhaps ballistic point 
contacts as 'nano contacts' might ultimately work, with particular few- 
A-thick semi-conducting nanospacers, and with very clean interfaces. 

The suggested experiment concerning the TMR would therefore be 
analogous to the very recent set-up of Garcia et ai, but with a 
small semi-conducting tip of the above-mentioned kind. 



5 Conclusions 

We have performed a systematic study of the CPP-GMR and CIP- 
GMR (i.e. with Current Perpendicular to the Planes or Current in 
the Planes, respectively), for systems consisting of two 'ideal leads' 
attached to two ferromagnetic metallic slabs, which are three mono- 
layers thick, separated by n s monolayers of non-magnetic (metallic 
or non-metallic) spacer material. A simple two-band model is used 
throughout, where in the metallic case one band mimics the s-bands, 
and the other one the d-bands which are spin-dependent in the ferro- 
magnetic metals, whereas for nonmetallic spacers the two bands de- 
scribe the valence band and the conduction band, respectively. The 
s-d-hybridization is taken into account, and is important for our re- 
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suits. For the metals, and also for the semi-conducting spacer, we use 
two specific parameter sets, corresponding in the metallic case very 
roughly to the situation of a Co/Cu/Co-trilayer system with Cu leads, 
whereas in the non-metallic case the two parameter sets of the spacer 
distinguish between situations corresponding to an 'indirect' and a 'di- 
rect' energy gap, respectively. 

The calculations have been performed by very accurate Green's 
function methods and applied to systems without impurities, with non- 
magnetic bulk impurities in the spacer, and with magnetic and non- 
magnetic interface impurities produced by mutual interdiffusion of the 
atoms near the interface. Within our formalism other situations could 
also have been treated, e.g. more complicated geometries as (i) the 'con- 



strictions' studied experimentally by Garcia el al, or (ii) quasi- 
one-dimensional problems, where the potentials are constant within 
a layer, but random from one layer to the next. The last-mentioned 



problem has already been extensively studied by J. Mathon, [171], and 



in two-dimensional systems by one of the present authors, ||42|| . Also 
the recent work of Sanvito et al. should be mentioned in this broader 
context. 

In all cases the influence of the strength of disorder leads to strong 
effects, sometimes to an enhanced GMR and sometimes to a reduction. 
However for different model parameters the outcome can be quite dif- 
ferent, which is not astonishing after all, since three conductances are 
involved in the changes, namely r| ++ ' ) , rj ++ ' ) , and r^ +_ ), and the GMR 
effect is - after all - a difference effect out of these quantities. 

So our calculation shows among other results that for the GMR one 
can hardly avoid extensive calculations for very realistic models of the 
specific systems considered. Only for the CPP case with impurities in 
the bulk of a metallic spacer, the GMR seems to be strongly reduced 
with increasing strength a T of the disorder, but only as long as a r 
remains < 1 (Fig. 8 and Fig. 9). In contrast, interface impurities may 
lead in the same case to a reduction of the GMR for one parameter set 
(Fig. 10), but to an enhancement for another set (Fig. 11). 

One of our main results is perhaps the suggested existence of a 
drastic increase of the GMR with increasing spacer thickness in case of 
ballistic conductance through a nonmetallic spacer, see section 4 above, 
where we also suggest an experimental realization along the lines of the 
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recent experiments of Garcia et a/., 
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Tables 



Model parameters (i), from overall band-structure 




Co 


Cu 




E(s) 


E(d,<7) 


E(s) 


E(d) 


T 


0.0 


-2.6 


0.0 


-2.6 


I 


0.0 


-2.1 


0.0 


-2.6 


Fermi energy: Ep = 3D — 2.8 




Hopping integrals: t s = 3D 


-l,t d 


= 3D - 0.2, t ad = 3D0.3 



Table 1: Model parameters (i); metallic spacer 



Model parameters (ii), from overall behaviour of the DOS 




Co 


Cu 




E(s) 


E(d,a) 


E(s) 


E(d) 


T 


0.0 


-1.0 


0.0 


-1.0 


1 


0.0 


-0.2 


0.0 


-1.0 


Fermi energy: Ep = 3D0.0 






Hopping 


integrals: t s = 3D 


-l,t d 


= 3D - 0.2, t ad = 3D1.0 



Table 2: Model parameters (ii); metallic spacer 





'Isolator 1' 


f a = -1.0,k=-0.2,t ad =0.3 


small gap : 0.025eV 


E s = -8.77 E d = -1.60 


intermediate gap : 0.17eV 


E s = -8.85 E d = -1.55 


large gap : 3AeV 


E s = -9.77 E d = -0.60 



Table 3: Parameter set (iii): Non-metallic spacer; 'indirect energy gap'; 
Ep = —2.8; the set of the first line was always used in the paper, with 
t s = —1, t d = —0.2, t sd = 0.3 throughout. 
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'Isolator 2' 


^ = -1,^ = 0^,^ = 0.3 


energy gap at fc[m] = 0.814 


l.OeV 


£ s = 2.50 £ d = -2.85 


energy gap at fc[m] = 1.888 


l.OeV 


E s = 0.00 E d = -2.80 


energy gap at fc[m] = 3.066 


l.OeV 


E s = -4.0 £ d = -2.80 



Table 4: Parameter set (iv): non-metallic spacer; 'direct energy gap'; 
Ep = —2.8; the arrow in Fig. 5 denotes the position of the minimal en- 
ergy difference E c (k) — E v (k); the set of the first line was always used in the 
paper; t s = — 1 and t s j = 0.3 throughout, as before, *) and = —0.2 for the 
hopping between two neighbouring metal atoms; whereas, if at least one of 
the two atoms were nonmetallic, this time = was assumed, so that the 
dispersion of the valence band in the r.h.s. of Figure 5 is solely caused by 
t s = — 1 and the on-site s-d-hybridization t st j_ = 0.3. 
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Figures 




Figure 1: Schematic plot of our systems for the CPP-geometry 




Figure 2: Schematic plot of our systems for the CIP-geometry 
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Indirect Energy Gap Direct Energy Gap 




k[in] k[ni] 



Figure 5: Band-structure of the non-metallic spacers (iii) and (iv); 
Ep = —2.8; on the left h.s. the band structure of model (iii) is shown, 
with an 'indirect energy gap' from the upper edge of the valence band 
at km = y/37i to the lower edge of the conduction band at k — 0, which 
is two times the energy distance to E F (see the inset); the r.h.s. is for 
model (iv) with a 'direct gap' represented by the arrow. 
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Figure 6: k-space calculation, CPP, system (i) 
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Figure 7: k-space calculation, CPP, system (ii) 
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Figure 8: x-space calculation, CPP, system (i), bulk impurities; the 
standard-deviation a r of Eq. ([]) characterizes the 'strength of disor- 
der'; 20 samples per point 
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Figure 9: x-space calculation, CPP, system (ii), bulk impurities; the 
standard deviation a r of Eq. @ characterizes the 'strength of disorder'; 
20 samples per point 
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Figure 10: x-space calculation, CPP; system (i); interface impurities; 
100 samples per point 
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Figure 11: x-space calculation; CPP; system (ii); interface impurities; 
100 samples per point 
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Figure 12: x-space calculation, CIP, system (i); no impurities 
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Figure 14: x-space calculation; 
'standard deviation' a r of Eq. ( 
der'; 100 samples per point 
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Figure 15: x-space calculation; CIP; system (ii); bulk impurities; the 
'standard deviation' o r of Eq. ([|) characterizes the 'strength of disor- 
der'; 100 samples per point 
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Figure 16: x-space calculation; CIP; system (i); interface impurities; 
100 samples per point 
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Figure 17: x-space calculation; CIP; system (ii); interface impurities; 
100 samples per point 
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Figure 18: CPP-GMR for non-metallic spacer (iv), i.e. with direct gap; 
NF\ = NF2 = 3; parameters corresponding to line 1 in Table 4; no 
impurities 
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Figure 19: CPP-GMR for non-metallic spacer (iii), i.e. with indirect 
gap; NFi = NF 2 = 3; parameters corresponding to line 1 in Table 3; 
no impurities 
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